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We study the effects of neutrino oscillation on supernova neutrino in the case of the 
inverted mass hierarchy (7713 >C mi < 7712) as well as the normal mass hierarchy (mi < mi -C 
7713). Numerical analysis using a realistic supernova and presupernova model allows us to 
discuss quantitatively a possibility to probe neutrino oscillation parameters. We show that 
£C) • information about the mass hierarchy can be obtained if #13 is rather large (sin 2 2#i3 > 10 - 3 ) 

^ ' and that #13 can be well probed by SuperKamiokande if the neutrino mass hierarchy is the 

inverted case. Errors due to the uncertainty of the original neutrino spectra and the Earth 
£ — ■ effect are also discussed. 
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Neutrino mixing and mass spectrum are the keys to probe new physics beyond 
the standard model of particle physics. Some of the neutrino oscillation parameters 
have been revealed dramatically by the observation of the atmospheric neutrino ^ 
and the solar neutrino 2 )" 6 ). Recently the first results of the KamLAND experiment 
have confirmed the Large Mixing Angle (LMA) solution of the solar neutrino problem 
7 ). An upper bound on #13 has also been obtained from CHOOZ experiment 8 ) and a 
lower bound is expected to obtained from single and double beta decay experiments 
9 ). But there still remain some ambiguities in the properties of neutrinos: the 
mass hierarchy, i.e., normal or inverted and the magnitude of #13. Current status is 
reviewed by many authors 10 ) " 13 ) . 

In such present situation, much attention have been paid to another neutrino 
source, supernova. This is a completely different system from sun, atmosphere, ac- 
celerator, and reactor in regard to neutrino energy, flavor of produced neutrinos, 
propagation length and so forth. Then neutrino emission from a supernova is ex- 
pected to give valuable information that can not be obtained from neutrinos from 
other sources. In fact, pioneering observations of neutrinos from SN1987A 14 )' 15 ) 
contributed significantly to our knowledge of the fundamental properties of neutri- 
nos 16 ) " 18 ) . Especially there have been many studies about the implication for the 
mass hierarchy from the observed neutrino events and the inverted hierarchy is dis- 
favored if sin 2 #13 is rather large (> 10~ 4 ) " 21 \ Here, normal and inverted mass 
hierarchies are the mass pattern 7713 3> 771,2 > m i and 771-2 > m i 3> 7773, respectively. 
In our notation Amfj = m 2 — mf, Z\m 2 2 and ^77i 2 3 ~ Am\^ are the mass squared 
differences which are related with the solutions of the solar and atmospheric neutrino 
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problems, respectively. There have also been studies to try to extract the unknown 
neutrino properties from future supernova 22 ) " 27 ) . 

In this paper, we calculate numerically the effects of neutrino oscillation on su- 
pernova neutrino, extending our previous study 27 ) where all analyses are performed 
with normal mass hierarchy, and investigate the possibility to identify the mass hi- 
erarchy and to probe the neutrino oscillation parameters #13 by the observation of 
the neutrinos from the next galactic supernova. We use the original neutrino spectra 
from supernova based on a realistic supernova model and the density profile of the 
progenitor star based on a realistic presupernova model. Since uncertainties of the 
original neutrino spectra are important in this analysis, we estimate the effect of 
them on our analysis. The Earth matter effect is also discussed, which have already 
been studied in the case of normal hierarchy in our previous papers 29 )> 30 ). 

This paper is organized as follows. In section II we summarize the properties 
of supernova neutrino briefly. The method of analysis is described and the results 
are shown in section III. We discuss some ambiguities in the basis of our study and 
summarize our results in section IV. 



§2. Supernova Neutrino 



Here we summarize the properties of supernova neutrino. For details, see, for 
example, a review by Suzuki 34 ). Almost all of the binding energy of the neutron 
star, 

GMl „ /M NS \ 2 /10km\ . , 

is radiated away as neutrinos. Here G, Mns and -Rns are the gravitational con- 
stant, the mass and radius of the neutron star, respectively. Due to the differ- 
ence of interaction strength, average energies are different between flavors. Al- 
though quantitative estimate of the difference is difficult, it is qualitatively true 
that (E Ue ) < (Ep e ) < (E Ux ). Here u x means v fl ,v T and their antineutrinos. These 
differences are essential in this paper. 

We use a realistic model of a collapse-driven supernova by the Lawrence Liver- 
more group 36 ) to calculate the neutrino luminosities and energy spectra, as we did 
in our previous paper 27 ). The average energy of each flavor is: 

(E Ue ) ~ 13MeV, (E Pe ) ~ 16MeV, {E U J ~ 23MeV. (2-2) 

Details of this original neutrino spectra are discussed by Totani et al. 37 ) These 
neutrinos, which are produced in the high dense region of the iron core, interact 
with matter before emerging from the supernova. Due to the nonzero masses and 
the mixing in vacuum among various neutrino flavors, flavor conversions can occur 
in supernova. When the mixing angle is small, these conversions occur mainly in the 
resonance layer, where the density is 

B . ( Am 2 \ /10MeV\ /0.5\ „ , , 

*- ~ x 10 g/cc [uv* ) ) ( n ) cos ^ m 
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Fig. 1. Schematic level crossing diagram for normal (left) and inverted (right) mass hierarchies. 
The circles shows the resonance points. 

where Am 2 is the mass squared difference, 6 is the mixing angle, E v is the neutrino 
energy, and Y e is the mean number of electrons per baryon. Since the supernova core 
is dense enough, there are two resonance points in supernova envelope. One that 
occurs at higher density is called H-resonance and another is called L-resonance. If 
the mass hierarchy is normal, both resonances occur in neutrino sector. On the other 
hand, if the mass hierarchy is inverted, H-resonance occurs in antineutrino sector 
and L-resonance occurs in neutrino sector. The schematic level crossing diagram for 
normal and inverted mass hierarchies are shown in Fig. 1. 

The dynamics of conversions including large mixing case is determined by the 
adiabaticity parameter 7, which depend on the mixing angle and the mass-squared 
difference between involved flavors: 

#13 and Am 13 at H — resonance, (2*4) 

612 and Am\ 2 at L — resonance. (2-5) 

When 7 > 1, the resonance is called 'adiabatic resonance' and the fluxes of the 
two involved mass eigenstate are completely exchanged. On the contrary, when 
7 <C 1, the resonance is called 'nonadiabatic resonance' and the conversion does not 
occur. The dynamics of the resonance in supernova is studied in detail by Dighe and 
Smirnov 22 ) . 

§3. Method and Results 

In this section we describe the method of analysis and show the results. 

3.1. Conversion Probabilities 

In the framework of three-flavor neutrino oscillation, the time evolution equation 
of the neutrino wave functions can be written as follows: 



H(t) 



(3-1) 
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(3-2) 



H(t) = U Am\ 2 /2E V 

v Amj 3 /2E U J 

where A(t) = v / 2G j pn e (t), Gf is Fermi constant, n e (t) is the electron number density, 
and Amjj is the mass squared differences. In case of antineutrino, the sign of A(t) 
changes. Here U is the Cabibbo-Kobayashi-Maskawa (CKM) matrix: 



U 



( C12C13 S12C13 S13 \ 

-S12C23 - C12S23S13 C12C23 - S12S23S13 S23C13 
\ S12S23 - C12C23S13 -C12S23 - S12C23S13 C23C13 J 



(3-3) 



where Sij = sin 9ij , Cij = cos6y for i,j = 1,2, 3(i < j). We have here put the CP 
phase equal to zero in the CKM matrix. 

By solving numerically these equations along the density profile of progenitor, 
which is calculated by Woosley and Weaver 39 ), we obtain conversion probabilities 
P(a — > /?), i.e., probability that v a at the center of the supernova becomes vp at the 
surface of the progenitor star. 

In our previous paper 27 ) , we assumed the normal mass hierarchy and took four 
models for neutrino oscillation parameters, the differences being the solution of the 
solar neutrino problem (LMA or SMA) and the magnitude of #13 . Here we take the 
following values: 

sin 2 26»i2 = 0.84, Am\ 2 = 7 x l(T 5 eV 2 , 
sin 2 26*23 = 1.0, Am 2 23 = 3 -2 x 10~ 3 eV 2 . (3-4) 

Values of #12 and Z\m 2 2 are taken from the global analysis of the solar neutrino 
observations and the KamLAND experiment 40 ) and correspond to the LMA solution 
of the solar neutrino problem while those of #23 and Z\m 2 3 are taken from the analysis 
of the atmospheric neutrino observation ^ . As to 6*13, we take two fiducial values 
sin 2 26*13 = 0.043, 10~ 6 as we did in our previous paper 27 ). Later we will discuss 
the case of the other values. Furthermore, we consider both normal and inverted 
hierarchy. Consequently, there are four models and we call them normal-LMA-L, 
normal-LMA-S, inverted-LMA-L and inverted-LMA-S. The last character (L or S) 
represents the magnitude of #13 (large or small). In our notation, Amfj = rnj — mf so 
that Am\ 3 > in the normal hierarchy case and Am\ 3 < in the inverted hierarchy 
case. Therefore, normal-LMA-L and inverted-LMA-L are different only in the sign 
of Aml 3 . 

We show in Fig. 2 demonstrations of conversion probabilities. The left figure 
is the time evolution of P(e — > e), probability that u e remains v e , and the right 
figure is that of P(e — > e), probability that v e remains v e . Four curves for the 
same model correspond to the neutrino of energy, 5 MeV, 10 MeV, 40 MeV and 
70 MeV, respectively. It can be seen that H-resonance and L-resonance occur at 
the O+Ne+Mg and He layer, respectively. In the neutrino sector, the conversion 
probabilities for the inverted hierarchy are the same as those for the model normal- 
LMA-S. This is because the H-resonance is completely nonadiabatic when #13 is very 
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Fig. 2. Time evolution of conversion probabilities P(e — ► e) and P(e — > e). Four lines of the same 
marking correspond to neutrino energy, 5MeV, lOMeV, 40MeV, and 70MeV, respectively. 



small, as in normal-LMA-S, and this is phenomenologically as if the H-resonance is 
absent as in the inverted hierarchy case. This logic also applies to why normal-LMA- 
L, normal-LMA-S and inverted-LMA-S are degenerate in the antineutrino sector, if 
we consider the H-resonance in the antineutrino sector. These degeneracies are the 
origins of the degeneracies that appear later in the event rates. 

3.2. Event Rates 

After obtaining the conversion probabilities, the neutrino fluxes at the Earth are 
calculated by multiplying the conversion probabilities by the original spectra and the 
distance factor l/4nd 2 . Here we take 10 kpc for the distance d between the Earth 
and the supernova. Further, by multiplying these fluxes by the cross sections of the 
detection interactions, the detector volume and the detector efficiency, we obtain the 
event rates at the detectors. Here we consider two detectors: SuperKamiokande (SK) 
and SNO. Properties of these detectors and cross sections used to calculate event 
rates are described in our previous paper 27 ). Unfortunate accident at SK lessened 
the detection efficiency at low energy (< 8 MeV) but this cause negligible effect in 
the subsequent analysis. 

Fig. 3-5 show the time-integrated energy spectra (left) and the time evolution 
of the number of neutrino events (right) at SK and SNO (u e charged current (CC) 
events and v e CC events), respectively. In Fig. 3, only v e p CC interaction is taken 
into account. Event numbers of each interaction are shown in Table I and II. In 
these tables, the contribution from neutronization burst phase is also shown. Here 
the neutronization burst phase means the period from 41msec to 48msec after the 
bounce. 

3.3. Distinction between Models 

In general neutrino oscillation makes the v e and v e spectra harder, since the 
original average energies of v e and v e are smaller than that of v x . In other words 
neutrino oscillation produces high energy v e and u e from v x . As a result, the high- 
energy events increase and the low-energy events decrease. The boundary between 
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Fig. 3. The time-integrated energy spectra (left) and the time evolution of the number of neutrino 
events (right) at SK. In these figures, only v e p CC interaction is taken into account. 
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Fig. 4. The time-integrated energy spectra (left) and the time evolution of the number of neutrino 
events (right) at SNO. In these figures, only v e d CC interaction is taken into account. 
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Fig. 5. The time-integrated energy spectra (left) and the time evolution of the number of neutrino 
events (right) at SNO. In these figures, only u e d CC interaction is taken into account. 



high energy and low energy is around 20 MeV. Note that how much these increase 
and decrease are depends on the adiabaticity parametes, and then the neutrino 
oscillation parameters, as can be seen in Fig. 3, 4 and 5. This feature can be used 
as a criterion of the magnitude of the neutrino oscillation effects. We define the 
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Table I. Number of events at SK 



hierarchy 


normal 




inverted 






model 


LMA-L 


T 1\ T A C 


T A \ A T 

LMA-L 


LMA-b 


no osc 


v e p 


9459 


9427 


12269 


9441 


8036 




186 


171 


171 


171 


132 




46 


46 


56 


46 


42 


v^e~ 


25 


26 


27 


26 


30 


Vy.e~ 


24 


23 


12 


23 


24 


v T e~ 


25 


26 


26 


26 


30 


v T e~ 


24 


23 


12 


23 


24 


Ov e 


297 


214 


297 


214 


31 


Ol7 e 


160 


158 


296 


159 


92 


total 


10245 


10114 


13084 


10129 


8441 


burst 


15.7 


16.7 


20.1 


16.7 


12.4 


Table II. Number of events (CC) at SNO 




hierarchy 


normal 




inverted 






model 


LMA-L 


LMA-S 


LMA-L 


LMA-S 


no osc 


v e d(CC) 


237 


185 


185 


185 


68 


i7 e d(CC) 


118 


117 


190 


117 


82 


total 


355 


302 


375 


302 


150 


burst 


0.6 


1.1 


1.1 


1.1 


2.1 



following ratios of high-energy to low-energy events at both detectors: 

_ number of events at 20 < E < 70MeV 
SK = number of events at 5 < E < 20MeV ' ^ ' ' 

_ number of events at 20 < E < 70MeV 
SNO = number of events at 5 < E < 20MeV ' ( ' ' 

Note that the energy range for the definitions of R^k and -Rsno are different from 
those in the previous paper 27 \ The plots of Rsk vs Rsno are shown in Fig. 6. In the 
left figure, we consider only v e d CC events at SNO for i?sNO assuming v e d CC event 
and v e d CC event can be distinguished completely. On the other hand, in the right 
figure we assume that v e d CC event and v e d CC event can not be distinguished at 
all and we sum v e d CC events and v e d CC events for -Rsno- The error bars represent 
the statistical errors. 

Note that v e flux and v e flux have essentially different information about the 
neutrino oscillation parameters. For example, inverted-LMA-L and inverted-LMA-S 
are distinguishable from v e events but are not from u e events. So it is more effective 
to distinguish between models if v e d CC events and v e d CC events at SNO can be 
distinguished perfectly. This can be clearly seen in Fig. 6. In the left figure it 
is easier to distinguish between normal-LMA-L and (normal-LMA-S and inverted- 
LMA-S) than in the right figure. But even in the left figure, it may be difficult 
to distinguish between normal-LMA-L and (normal-LMA-S and inverted-LMA-S) 
considering some ambiguities discussed in the following sections. 
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Fig. 6. Plots of Rsk vs -Rsno- In the left figure, only v e d CC events at SNO are considered for 
J?sno- In the right figure we sum v e d CC events and v e d CC events to calculate -Rsno- The 
errorbars show statistical errors only. 



3.4. Uncertainties in the original spectra 

One of the crucial ingredients in this study is the flavor-dependences of the 
original spectra. Smaller differences will make it more difficult to distinguish between 
various models. There have been many numerical simulations of core collapse but 
the predicted luminosities and temperatures of neutrinos are different from group by 
group. Numerical model by the Livermore group 36 ), adopted in this paper, has the 
great advantage that it covers the full evolution of supernova: from the core collapse 
over the explosion to the cooling phase of the protoneutron star, although it involves 
rather traditional treatments of neutrinos. 

There are simulations with more sophisticated treatments of neutrinos but they 
do not obtain explosions and neutrino fluxes for only less than 1 sec are available 
now 41 )' 42 ). These simulations predict less average-energy differences between flavors 
compared to those of the Livermore group: E x /E e ~ 1.4 41 \ 1.7 42 \ E x /E s ~ l.l 41 ), 
1.3 42 ), where E e , E s and E x are the average energies of v e , v e and u x , respectively. To 
see the effect of the uncertainties in the temperature differences on our analysis, we 
perform similar analyses varying v x temperature as stated below and fixing v e and v e 
temperatures. This is relevant because various simulations agree well with respect 
to v e and 9 e temperatures. The original spectra can be fitted by the "pinched" 
Fermi-Dirac distribution, 

f(E,T a ,r la )<X eE/Ta E \ ria+v (3-7) 

where the temperature T a and the pinching parameter rj a for each flavor (a = e,e 
and x) are 

(T a ,rj a ) = (2.8MeV,2.6) for u e , 
(4.0MeV,2.5) for v e , 

(7.0MeV,0) forzv (3-8) 
We consider the following transformation of the distribution function of v x which 
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show statistical errors only. 



has temperature T x , 

f(E,T x , Vx ) - (ffl HE^XiV*) (3-9) 

and we regard it as the distribution function of v x which has temperature T' x . Note 
that this transformation conserves the total energy of v x . Based on this spectrum 
with various values of T x , we obtain Rsk and -Rsno m the same way as in the 
previous subsections. Fig. 7 shows T x /T e and T x /T s dependencies of -Rsno and 
Rsk, respectively. In the left figure, only v e events are considered to calculate -Rsno- 
Predicted values by the Livermore group simulations are T x /T e = 1.9 and T x /T s = 

1.4, respectively. It should be noted that T x /T e and T x /T s are not independent 
parameters here since we fix T e and T s . 

As can be seen, smaller values of these ratios result in smaller separation of event 
ratio between different models. In the left figure the errorbars overlap each other if 
the temperature ratio is less than 1.5. On the other hand, it seems that inverted- 
LMA-L and the others can be discriminated independent of the temperature ratio. 
The reason why Rsk is significantly different between inverted-LMA-L and the others 
is that the pinching parameter is different between v e and v x . By the same reason, 
-Rsno f° r normal-LMA-L is not the same those for the others even if T x /T e = 1. 
Although the relevance of the transformation (3-9) is difficult to argue when T x /T s 
is nearly unity, the possibility of the discrimination between inverted-LMA-L and 
the others will be robust against the change of v x temperature of about 10%. 

3.5. Dependence on 0i 3 

So far we have considered only two extreme cases: H-resonance is perfectly 
adiabatic and nonadiabatic. It is interesting to investigate the intermediate cases. 
In Fig. 8 we show 6> 13 dependence of i?sK and i?sNO Only v e events are taken 
into account to calculate -Rsno Note that Rsk and -Rsno vary only in the case of 
inverted and normal hierarchy, respectively, as will be expected from Fig. 6. In the 
case of normal hierarchy it will be difficult to determine the value of # 13 due to large 
statistical errors but will be possible to say whether it is very large or very small. 
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On the other hand, in the case of inverted hierarchy the overlap of the errorbars are 
small even in the intermediate cases. If #13 is rather large (sin 2 2^13 > 10~ 3 ) the 
mass hierarchy will be identified. 

3.6. Earth Effects 

Depending on the detector position, neutrino go through the Earth and the 
matter effect inside the Earth can change the neutrino spectra again. For the matter 
density of the Earth, neutrino oscillation between only two light neutrinos is involved. 
Since the oscillation length in the Earth, which depends on the neutrino energy, is 
the same order as the Earth radius for the neutrino parameter adapted here, the 
Earth matter effect appear as a distosion in the spectra 29 ) " 33 ). 

The numerical process to calculate the Earth effect is the same as described 
above except that the density profile of the Earth is needed. We use a realistic 
density profile 43 ) . In Fig. 9 neutrino spectra in the case that the nadir angle is 30° 
and 70° are shown with that in absence of the Earth effects. The left figure is the 
u e spectra at SNO and the right figure is the v e spectra at SK. The Earth effect is 
absent in v e and v e in case of normal-LMA-L and inverted-LMA-L, respectively. This 
is because in these cases H-resonance is perfectly adiabatic and low-energy neutrinos, 
which were originally v e in the supernova core, are converted to the heaviest neutrino, 
which is not involved in the matter oscillation in the Earth. Thus the detection of 
the Earth effect will be helpful to distinguish models, especially normal-LMA-L from 
normal-LMA-S and inverted-LMA-S. 

The form of the distorted spectra depends on the nadir angle of the neutrino 
path. The nadir angle 30° and 70° correspond to one of the path with which neutrino 
pass through core and mantle, and only mantle, respectively. The detectability of the 
Earth effect was discussed by us 30 ) and it is shown that complementary observation 
by SK, SNO and Large Volume Detector (LVD) is effective for its detection. 
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§4. Discussion and Conclusion 

There are some ambiguities besides those descussed above. The first is the 
direction of the supernova. If the supernova can be observed optically, the direction 
can be known with enough accuracy. But if the supernova is at the Galactic Center, it 
might be hidden by the large amount of gas and could not be seen optically. Pointing 
by the electron scattering events of the supernova neutrino is studied by several 
authors 44 )' 45 ) and the accuracy is expected to be ~ 7°. More detailed analyses of 
the Earth effects have been studied considering the locations of the detectors and 
the direction of the supernova 30 )> 32 ). 

Another is the mass of the progenitor star. It affects the mass of the iron core, 
which affects the neutrino spectra 42 )' 46 )' 47 ). Study including the mass uncertainty 
is now in progress but the preliminary result is that the mass uncertainty is not 
important in our analysis 48 ). 

Recently effects of shock propagation on neutrino oscillation in supernova have 
been studied 49 )' 50 )' 24 )' 51 ) and it was shown that some characteristic signatures may 
emerge as the shock propagates through the regions where matter-enhanced neutrino 
flavor conversion occurs. As we show 50 ), shock propagation effect will be safely 
removed by taking only early-phase events into account. 

We studied the effects of neutrino oscillation on supernova neutrino in the case 
of the inverted mass hierarchy as well as the normal mass hierarchy. Numerical 
analysis using a realistic supernova and presupernova model allowed us to discuss 
quantitatively a possibility to probe neutrino oscillation parameters. We showed that 
degeneracy exists only between normal-LMA-S and inverted-LMA-S if the Earth 
effect is taken into account and that 6*13 can be well probed by SK if the neutrino 
mass hierarchy is inverted case. Errors due to the uncertainty of the original neutrino 
spectra are also estimated. 
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